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Abstract 

Groundwater age has become a fundamental concept in groundwater hydrology, 
but ages originating from isotopic analyses are still identified with a lack of 
clarity and using models that occasionally are unrealistic. If the effect of advec- 
tion and dispersion on water ages has already been extensively identified, very 
few studies address the reliability of using radiometric ages as derived from 
isotopic data to estimate aquifer properties such as average velocities. Using 
simple one-dimensional and two-dimensional analytical solutions for single-site 
and two-sites mobile-immobile systems, we compare the radiometric ages to the 
mean ages (or residence times) as deduced from a direct, physically-based simu- 
lation approach (using the mean age equation), and show that the competition 
between isotope decay rate and dispersion coefficient can generate important 
discrepancies between the two types of ages. A correction for the average ap- 
parent velocity originating from apparent isotopic ages is additionally provided. 
The particular case of the Tritium age dating method is also addressed, and a 
numerical example is finally given for illustrating the analysis considering a more 
complex and heterogeneous aquifer system. Our results suggest that age defi- 
nitions based on the radioactivity of isotopes may not be representative for the 
mean age of the sample or for the groundwater velocity at given locations, and 
may not always be suitable for constraining the calibration of hydrogeological 
models. 

Keywords: Age dating. Apparent age. Mean age. Tritium age. Apparent 
velocity 



1. Introduction 

Radiometric age dating has been popular for decades in groundwater hydrol- 
ogy and for more than a century in other disciplines such as geology, archeology. 
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and so on. It is one of the recent advances in hydrogeology to clarify the defini- 
tions of the kinematic, mean, and radiometric ages of groundwater. In general, 
there is an increasing consensus that the apparent radiometric age, which is 
the most measurable quantity, cannot be considered to represent the mean age, 
which may be the most representative for the average residence time of water 



molecules in the subsurface (Bethke and Johnson (2002), Bethke and John 



son 



(20081, Sanford (2011)). Discrepancies between ages obtained by means 



of different approaches have been frequently reported, and in the majority of 
the cases purely advective ages are compared to ages definitions that include 



advection and dispersion processes (Castro and Goblet (20051, LaBolle et al. 



(2006), Bethke and Johnson (2008)). Concentration-based ages are simply ap- 



parent ages that do not a priori relate to the mean age of a water sample ( San 
ford| ( |2011[ )). Radiometric ages are classically derived from the interpretation of 



chemical/isotopic data that ignores transport characteristics although isotope 
transport undergoes dispersion/diffusion, decay and to some extend reaction 
processes in the subsurface. The kinematic age of groundwater is defined as 
the time for which a groundwater molecule may travel along a one-dimensional 
stream line from the time of recharge (birth). Thus, this kinematic age is calcu- 
lated based only on the advective component of the migration of water molecules 
(piston flow) without considering the mixing processes among water molecules 
which enter the system at different times in the past (different birthdates). Since 
the kinematic age does not consider mixing, it cannot be directly translated from 
any chemical signature in the real world where various types of mixing always 
exist. Irrespective of the importance of the kinematic age it should be noted that 
it is a function of only the groundwater velocity which is one of the most fun- 
damental physical parameters to describe the groundwater flow system. When 
the ages are known at two different locations, hydrogeologists may hope to be 
able to calculate the average velocity of groundwater dividing the distance by 
the difference in age along this distance. 

A radiometric age is a measurable quantity for a given water sample, based on 
the change in the activity of decaying solutes. When the initial activity at birth 
and the rate of the decay (change) are known, the radiometric age can be calcu- 
lated from the current activity at given locations. This assumes that the water 
sample under consideration is in a closed system, i.e. a system that does not 
exchange any material with its surrounding. Radiometric age thus differs from 
the mean age of groundwater, which is defined as the first temporal moment 
of the age probability density function for given samples, describing how age 
distributes as a result of advection and mixing processes. Age mass is a useful 
extensive (additive) property, defined as the product of age and mass for a given 
mass (e.g. see |Goode| (|1996[), |Varni and Carrera] (|1998[), |Ginn| (|1999|). Math- 



ematically an intensive property or the age probability density, defined as the 
age mass divided by the total mass, is a popular concept to derive the equation 
that describes the mixture of different ages of groundwater. 



Bethke and Johnson ( 2008j ) clarify the difference between the sample apparent 



(radiometric) age and the mean age of groundwater, the mathematical descrip- 
tion of which has been well established by Goode ( 1996 ) , iVarni and Carrera 
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(1998), Ginn (19991, Cornaton and Perrochet (20061, and Ginn et al. (20091 



among others. Although it seems that there are many different versions of 
mathematical equations to describe the mean age, age probability density, and 
age mass, conceptual difference is minor and consists in the assumed steady or 
transient flow conditions or different dependent variables as an integral moment 
or probability density (Ginn (2007), Cornaton and Perrochet (2007)). 
The recent finding that the radiometric age may not be representative for the 
mean residence time or the velocity of groundwater can imply that the only way 
to analyze the mean residence time is to consider applications of the mean age 
equation. This may lead to questioning the relevance of the radiometric ages 
frequently estimated by isotope hydrologists and for what purpose the radio- 
metric age is measured. Such an age may not be the most relevant data to be 
used for calibrating a groundwater flow field in a numerical model (for example 
by solving inverse problems constrained by radiometric age data and traditional 
hydrological data). As indicated by Bethke and Johnson (2008), time is ripe for 
developing a new paradigm that would account for the recent findings. 
In this paper, we first review mathematical equations to describe the transport 
of the kinematic age, mean age, and decaying solute and illustrate the differ- 
ences among them. The relationships between a measurable quantity and more 
intrinsic quantities are then analyzed. A discussion on what can be interpreted 
from measurable quantities and under which conditions is finally made. 



2. Mathematical Description of Kinematic, Radiometric and Mean 
Ages 

This Section is dedicated to the definition of the kinematic, radiometric and 
mean ages, and to the comparison between each age property. One-dimensional 
and two-dimensional analyses are made in order to reveal the bias induced by 
the use of apparent ages for calculating apparent flow velocities. 

2.1. One- Dimensional Analysis 

Considering a simple one-dimensional flow system where groundwater recharges 
at the origin [x = 0) and the flow system is homogeneous in a steady state 
regime {v{x,t) = w), the kinematic age {Ak) is straightforward to calculate 
from its definition: 

Ak^- (1) 

V 

Steady-state transport of a decaying solute can be described using a simple 
advection-dispersion-reaction equation as: 

ax ax'^ 

where D is the dispersion coefficient and A is the constant for the first order 
decay to satisfy C{t) = C{t = 0)e^'*'*. Since the boundary concentration at 
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origin is assumed to be known, C{x = 0) = Co, the solution of Eq. ^ is 

2D 

where 



CM = Coexpf ""'' ^' 'l (3) 



and where use has also been made of the condition — D'^H^ = 0. The 

x—oo 

apparent age {Aa) is a direct translation of the initial concentration and the 
rate of change in the first order kinetic decay and is classically given as 

It is obvious from Eq. ([5| that the apparent age is influenced by the mixing 
{D) and the solute decay rate (A) as well as by the flow velocity [v) and the 
location {x). Eqs. ([T]) and ^ also indicates that apparent age Aa can only re- 
late to kinematic age for the asymptotic situation, i.e. Aa = Ak when ^ — > oo. 



The equation describing mean age [Am) is well established (e.g. see Goode 



v^^^D—^ + l = Q (6) 



(1996), Varni and Carrera (1998)) and a simplifled version for the given one 
dimensional domain can be written as 

dAjyi(x) ^ ^d^Ajyi[x 
dx dx 

When the birth for groundwater is assumed at the inlet and a boundary condi- 
tion for the mean age is given as Am{x = 0) = 0, the mean groundwater age is 
same as the kinematic age: 

Am{x) = ^ (7) 
for which the condition — D ^^q}^^ — has also been used. Equations ( 1 1 

a:— oo _ 
imply that if the flow system is one-dimensional or if there is only one 

velocity component, the mean age which is same as the kinematic age is only a 

function of a flow velocity and of the location in the flow system, and thus it is 

straightforward to translate the mean age to describe the flow velocity. In other 

words, when one knows the groundwater mean ages {A„i{xi) and Am{x2)) at two 

different locations {xi and X2) in a one-dimensional flow system, the distance 

between two locations divided by the difference in mean age can provide an 

average velocity between two points: 

(8) 



Am{x2) - A^{xi) 



When the apparent radiometric ages {Aa{xi) and Aa{x2)) are known, the aver- 
age velocity between two locations can also be calculated from ^ such that 

V = Va (9) 
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where the apparent velocity Va is defined as [x2 — Xi]/[Aa{x2) — Aa{xi)], by 
analogy with Eq. ([8|. Eq. ^ indicates that the groundwater velocity can be 
calculated using the apparent velocity when the decaying property of the solute 
(A) is known and the mixing processes (D) are clearly defined. It is also implied 
from ([9]) that the apparent velocity generally overestimates the actual velocity 
and that the difference between them becomes greater when Va is smaller and 
the product XD is larger. When the dispersion coefficient is defined a,s D = 
Ddiff + oiiv, Eq. ([9| can be re-arranged to give 

^ ^ Va{'^-^Dd^ff/vl) 
1 + Xai/va 

where Ddi ff is the effective diffusion coefficient in the porous medium and ai is 
the coefficient of longitudinal dispersivity. Defining the dimensionless variab le 



Ua = jj) (defined as Jenkins number in Mouchet and Deleersnijder (2008)), 
Eq. ([9]) takes the form: 

The variable IIci is similar to the Peclet number (Pe) if one considers a charac- 
teristic distance L = ^ as being the distance that tracer molecules can travel 

if they move at the average velocity Va during time A^^ (Pe = = ^ = Ha). 

Finally, the mean age of a radioactive tracer (A^.x) can be obtained by 
calculating the first two temporal moments of the corresponding age concentra- 
tion density function (pdf). The age concentration pdf g\{x, r) of a radioactive 
tracer (where r denotes the variable age) is solution of the advection-dispersion- 
reaction equation 

dgx{x,T) dgx{x,T) d^gxjx.T) 

with the boundary conditions g\{0,T) — S{t) (zero age flux concentration) and 



D 



0: 



^) = Tl^l^iT^ (- 4Dr^ - =exp(-Ar)5o(a:,r) (13) 

For instance, the zero-order temporal moment of g\{x,T) gives the mass of 
radioactive tracer in the system: 

r ^ r , (xv{l-(i)\ 
mo^x{x) = J gx[x, Tjdr = exp I — 1 (14) 

in which /3 has been defined in Eq. Q. The first-order temporal moment of 
g\{x, t) gives the mean age concentration of the radioactive tracer: 

mi,x(x) = / Tgx{x,T)dT ^ —mo^xix) (15) 
Jo 
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The mean age of the radioactive tracer is then obtained as 



f3 



(16) 



Note that mQ_x[x) satisfies the tracer mass equation ([2| using the boundary 



conditions too.a(O) = 1 (unit mass input) and — -D ^^q^^^^ 
mi, A (a;) satisfies the tracer mean age concentration equation 



— 0, and that 



dmi^xix) , d'^mi^xix) 



dx 



D- 



dx"^ 



Ami,A(x) + mo,A(2:) = 



(17) 



with the boundary conditions mi,A(0) = (zero mean age concentration) and 



D 



dx 



= 0. 



Dimensionless variables can be used to describe the relations between the 
apparent age, mean tracer age and mean water age, and between the apparent 
velocity and the actual velocity. The ratio of the apparent and mean ages 
ta = Aa/Am and the ratio of the apparent and mean velocities — Va/v 
can be defined using the dimensionless variables IT = v^/\D, Ha ~ 
^a.difsj^ vl/XDdiff,a.iid Ua^isvj^ Va/Xai. The following relations are derived 
from III, ([§, and ([lO]) (see Fig. [I]): 




1 - l/Ua,diff 



(18a) 



(18b) 



A first fundamental consideration arising from the analysis of Eqs. ( |16[ ), ([7| 



is that the inequality A^^x ^ Aa < A^ always holds, as already 



and (18a 
shown by 

this inequality is valid for more general contexts. Another interesting aspect is 



Deleersnijder et al. ( 2001 ) . Delhez et al. ( 2003 ) also demonstrated that 



that the age and velocity ratios in Eqs. (18a) and (18b I are independent on the 
x-coordinate. Eq. (18a) and Fig. [T^ indicate that the apparent age generally 
underestimates the mean age and that the ratio between two ages becomes larger 
as the decay and dispersion dominate the advective transport and vice versa. 



Eq. (18a) also indicates that the age ratio does not depend on the spatial 
coordinate x. When the parameter 11 decreases (i.e. XD increases), then the 
bias increases, and vice-versa. Thus apparent ages are reliable only if 11 is large 



enough, or consequently only if the product Al? is small enough. In Eq. ( 18b I 
and Fig. [TJd, the apparent velocity is similar to the actual velocity when the 
transport regime is advection-dominated (H > 1), and the actual velocity can 
be calculated from the apparent velocity when Ha > 1. When the dispersion 
dominates the transport {Ha < 1), it may not be possible to estimate the actual 
velocity from the apparent velocity. 
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Figure 1: The ratios of (a) apparent and mean ages, and (b) apparent and mean 
velocities as functions of the dimensionless variables 11 — /\D and Ha — v^/\D. 



2.2. Coupled One- Dimensional Solutions for Two- Dimensional Analysis of Aquifer- 
Aquitard Systems 

A simple two-dimensional domain can be considered when mobile and im- 
mobile (or a lot less mobile) waters co-exist in a groundwater system. Such 
two-dimensional analyses have been widely adapted for a fracture-matrix sys- 



tem (Cornaton et al. (2008)) and also in an aquifer-aquitard system (Bethke and 



Johnson] ( |2002[ |2008|), |Ginn et al.| (|2009|). In these flow systems, the flow in a 
mobile region is essentially one-dimensional and the velocity is large compared 
to that in an immobile area so that it is safe to assume that the velocity in 
the immobile region is negligible. When an aquitard (or a matrix block) has a 
dimension of 2L perpendicular to the flow direction x between two aquifers of a 
dimension 21 (or two fractures of an aperture value of 21), transport of the mean 
age or decaying solute can be described using an advection-dispersion-reaction 
equation in the aquifer and the diffusion-reaction equation in the aquitard (see 
Figure [2]). 

The kinematic age in the mobile aquifer (Ak.m) is, by definition (see Eq. (|1 
a direct translation of the velocity in the mobile region (vx)- 



Afi^rn. 



(19) 



Since the velocity in the immobile region is assumed to be negligible, the kine- 
matic age in the immobile aquitard {Ak^im) is infinite for all x and z, where z 
represents a direction perpendicular to the flow direction along the aquifer. 
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Figure 2: Illustration of tlie aquifer-aquitard system used to derive two-dimensional 
analytical solutions. 



Steady-state transport of a decaying solute in the aquifer is described by the 
following one-dimensional advection-dispersion-reaction equation such that 



dx 



Ur — 7, ^^m\^) + 



dx"^ 



dz 



(20) 



where and represent the porosities in the mobile and immobile regions, 
respectively. Steady-state transport in the aquitard can also be described by a 
one-dimensional diffusion-reaction equation: 



D 



d'^Ci„i{z\x) 



\CMx)=Q 



(21) 



When the concentrations in the interface between mobile and immobile regions 
are assumed to be identical {Cim{z\x) = C„i{x)) and the mass flux in the 
middle of the aquitard is negligible (_£)^f^^ii^iMi£) 
of Eq. (21 ) is given as 



— 0), then the solution 



■iz-L) 



-iz-L) 



(22) 



Using Eq. (22), Eq. (20) can be re- written as 

dCmix) , „ d'^Cmix) 



dx 



dx"^ 



iX + j)C^{x)=0 



(23) 
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where 



(24) 



The solution of Eq. ( 23 ) is similar to Eq. (|3| where A is replaced by (A + 7) 
Cm{x) = CTO(0)exp 



XVx 

2D, 



4(A + 7)i?, 



(25) 



and the apparent age in the aquifer {Aa^m) is also similar to Eq. ([s]) such that 



4(A + 7)i?, 



2AD, 



(26) 



The apparent age in the aquitard (^a.irn) can also be calculated using (22 1 
and pll): 



Aa,m{x) + ^ In 



(z-L) 



(27)- 

The results (261 and (27 1 indicate that the apparent age in the aquifer is de- 



termined by the decay rate of the solute and the mixing processes in both the 
aquifer and aquitard, and that the apparent age in the aquitard is same as that 
in the aquifer at the interface between the two domains and has a maximum 
value at the center of the aquitard when z — L. The maximum of the apparent 
age is found at z = L: 



Ta&-x[AaAm{x., z)] = Aa,im{x, Z = L) = Aa-m + ^ hi 



(28) 



The maximum age increases as the aquitard becomes thicker (larger L), the 
decaying rate of the solute is higher (larger A), and the effective diffusion coef- 
ficient {Dz) is smaller. Fig. |3] illustrates the distributions of the concentration 
of a decaying solute and the apparent age in such a two-dimensional system 
when I = 5m, L = 10m, — lOm/yr, A — ln2/100yr, — otiVx = lOOm^/yr, 
Dz = 0.004m^/yr, and rim = f^i™. = 0.1. 

Mean groundwater age in the aquifer {Am.m) and in the aquitard (^m,im) 
can be described using equations similar to ([6|, ( pO| and (21 ): 



dAjyi m (X) 

dx 



(fAr. 



,1 



dz 



+ 1 = (29) 
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Figure 3: Illustration of the distributions of (a) relative concentration of a decaying 
solute and (b) apparent age in a two-dimensional aquifer-aquitard system. 



1 = 



By applying the following boundary conditions: 







Am,m{x — 0) — 



The solutions of ( 29 ) and ( 30 1 are given as 



(30) 

(31a) 

(31b) 
(31c) 

(32a) 



-^m,imiz\x) — z{2L z) + „i(x) 



(32b) 



Eq. ( 32a ) corresponds to the result of Bethke and Johnson ( 2002 ) , the so-called 



groundwater age paradox, which is actually a simple direct consequence of the 
hypotheses of steady-state conditions, constant porosity and no diffusion in 
the direction of aquifer flow in the aquitard ( |Ginn et al. (2009)). Clearly, 



the mean age in the aquifer is a function of the rati o of the water volumes 
of the two systems ( f ) to the kinematic age. Eq. (32b) indicates that the 
age in the aquitard is a travel time by pure diffusion irom the aquifer to the 
aquitard midpoint. It is parabolic with its maximum located at the midpoint 
(max[Am, im] = Am,m + i^/2I?z, whcu z = L). Fig. [4] illustrates the distribu- 
tions of the groundwater mean age for the same groundwater flow system in 
Fig. [3] Fig. [5] illustrates the distributions of the kinematic, apparent, and mean 
ages along the flow direction in the mobile region and the apparent and mean 
ages across the immobile region (at x = 1000m) of the two-dimensional system 
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shown in Figures [3] and [4] It is clear from Fig. [5] that the apparent age over- 
estimates the kinematic age but underestimates the groundwater mean age in 
the aquifer. Groundwater mean age in the aquitard is significantly higher than 
the apparent age although both ages are influenced by the mixing processes in 
the aquitard. This is because the apparent age is influenced by the decay rate 
of the solute. From Eq. ( 26 1 , although the relationship between apparent age 
and mean age is more complex than in the one-dimensional system, it is unique 
and thus the apparent age can be used to estimate the mean age for a given 
distance from the source, groundwater velocity, the rate of decay for the solute, 
and the thickness and diffusive properties of the immobile region surrounding 
the mobile region. 




Figure 4: Illustration of the groundwater mean age distribution in a two-dimensional 
aquifer-aquitard system. 



2.3. Two-Dimensional Analysis 

Effects of transverse mixing can be analyzed by considering a two-dimensional 
system where the flow velocity can be approximated by a one-dimensional uni- 
form velocity (vx) but where recharge is limited at the origin [x — 0, y = 0). 
The kinematic age is approximated along the streamline from the source (y — 0) 
using the location and the one-dimensional velocity: 

Ak{x,y = 0)^— (33) 
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Figure 5: Illustration of the kinematic, apparent, and groundwater mean age distribu- 
tions (a) along the flow direction in the mobile region and (b) across the immobile 
region (x = 1000m) of the two-dimensional system illustrated in Figures (ll and (|4|. 



The steady-state transport of a decaying solute is described by 
dC{x,y) 



dx 



V • D • VC(a;,y) - XC{x,y) = 



(34) 



where the two-dimensional operator is defined as V = {d/dx,d/dy). When the 
dispersion coefficient D = {D^x, D^y] Dyx, Dyy) is such that Dxx — onVx+Ddiff-, 
Dx 



Dyx = 0, and Dyy = atVx + Ddiff, Eq. (34 1 simphfies in 



dC{x,y) 
dx 



Dx 



d^C{x,y) 
dx"^ 



D„ 



d^C{x,y) 
' dy'^ 



XC{x,y) = 



(35) 



For the given boundary condition of continuous mass injection with the rate of 
QCq (MT~^) at the origin, where Q is the inflow rate and Cq the flux concen- 
tration, the solution of Eq. (35 1 is given by 

QCo ' ' ' 



C(x,2/) 



2iniyJ DxxD. 



: exp 



vv 



2Dx 



ADx 



+ A 



Dx 



y_ 

D: 



(36) 

where n is the porosity and Kq is the modified Bessel function of second kind 
and zero order. When an asymptotic form of Kq{p) ( » ^yT r/2pe~P for large p) 
is introduced in Eq. ([36| (e.g. see Wilson and Miller (1978)), then it reads: 



C (X, y) — — ^^^^^^^^= exp 
n^^TiDxxDyyr/B 



X — r 
B 



where 



B = 



2D XX 



and 



2B\ 



9 Dxx o 

x- + -jz—y^ 

^yy 



(37) 



(38) 
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If the apparent age is defined from the ratio in concentrations of decaying and 
conservative solutes (A = 0), then it is given as 



C 



-In , 



1 /I 



A V4 



In K 



r 
B 



with 



K = 1 



2B\ 



(39a) 



(39b) 



(40) 



An asymptotic solution for A„i{x,y) is obtained when solving the mean age 
equation 

dAm{x,y) , ^ d'^Am{x,y) d'^Am{x,y) 

Vx ^ V iJxx TT^ \- Uyy TT^ h 1 = (41) 



dx 



dx^ 



9y2 



Using the zero mean age boundary condition Am{x — 0,y ^ 0) =0, this solution 
is: 

1 



Am{x,y) 




Amix,y = 0) 



(42a) 



(42b) 



Eqs. ( 39a ) and ( 42a I show a dependency of age with transverse dispersion Dy 



The ratio of apparent age and mean age is more complex than in the one- 
dimensional system but it still can be computed and potentially used to assess 



the dispersive properties of the aquifer. Eqs. (39b) and (42b) indicate that the 



apparent and mean ages are not influenced by the transverse dispersion along 
the streamline from the source. The apparent and mean ages, compared to the 
kinematic age, yield the following relation: 



TA^Ax.y^ 0) = 



Ar. 



= 1 



D 



XX 

XVx 



rA^^{x,y^ 0) 



Aa 

Ak 



4Ax 



In K 



2AD, 



(43a) 



(43b) 



In Eq. ( 43a ) , it is clear that ratio between the kinematic and mean ages becomes 
unity for large x and Vx-, and small Dxx- When Dxx ~ onVx, '^Am^ ~ \ -\- ai/x 
and thus the mean and kinematic ages becomes similar as the location is farther 



from the source compared to the longitudinal dispersivity. In Eq. ( 43b ) , it is also 

.2 

clear that ta 



(^ 



1) for large x, in which case the one-dimensional 



solution (18a) is recovered. As a result, it can be concluded that the effect of 



dimensionality on age bias is minor if the velocity field is uniform. 
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3. Implications for Tritium Age Dating 



Transient transport of isotopic signatures has been widely used to measure 
the age of groundwater. One of the most popular radioisotopes in dating mod- 
ern groundwater is the tritium which was produced by the atmospheric testing 
of thermonuclear bombs between 1951 and 1980 (activity peaked at 1963, e.g. 
Maloszewski and ZubCT| ( |1982[ ), rZuber et al.| ( |2005[ ), |Zuber et al.| ( |2001[ ), |Mor 



see 



genstern et al. (2010)). There can be various different approaches to dating 



groundwater using tritium as explained by Clark and Fritz ( 1997 ) and they 



may be categorized as: (i) identifying the spatial location of concentration (or 
activity) peak, (ii) calculating the time for decay from a known input level to 
the measured level, {in) analyzing concentration time series at a specific point. 
Clark and Fritzj (1997) indicated that the tritium age dating could be most 
appropriate when the influence of the mixing and dispersion is minimal for iso- 
tope activities. It is interesting to note that the spatial location of the tritium 
activity peak in the subsurface is influenced only by the advection and thus it 
is related only to the kinematic age, while the attenuation and decay can be 
strongly influenced by the mixing or other dilution processes. Both of these 
measures, however, could be misleading if they are used to analyze the mean 
residence time or mean age of the sample. 



Transport of tritium in a one-dimensional domain can be described by the 
following advection-dispersion-reaction equation: 



dC(x,t) dC{x,t) 

= —V- 



dt 



dx 



D 



d^C{x,t) 
dx^ 



XC{x,t) 



(44) 



For the simplistic analysis, a steady flow field and an initial Gaussian slug input 
of tritium is assumed at the origin: 



C(x,t = 0) 



1 



: CXp 



(45) 



The variable to denotes the initial time at which the concentration defined in 



Eq. ( 45 ) is taken as the initial concentration. The solution of Eq. ( 44 1 is given 

by 

1 / (x - vt)^ 



Cix,t) 



: exp 



-Ht + to) 



(46) 



y/4TrD{t + to) V 4D(t + to) 
In Eq. (46), the position of the peak is located aX x ~ vt with the peak value 

g-A(t+to) 



being 



C{x = vt,t) = 



^AirDit + to) 



(47) 



If the location of the peak (x — Xpeak) is known at given time t = tp^aki 
groundwater velocity is calculated as v — Xpeak/tpeak and the kinematic age of 
groundwater at x — Xpeak is given by Ak{x = Xpeak) — tpeak, which is same as 
the mean age of groundwater in a one-dimensional flow system [A^ = A.„i). It 
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is noted again that this is based only on the fact that the peak migrates with 
the rate of groundwater velocity and is not related to the decay of the solute. 

If the age is calculated from the change in peak concentration over the time 
period tpeak (also known as time lag method) , the apparent age of groundwater 
is given by 

Because the apparent age calculation is based on the rate of decrease in con- 
centration values, additional dilution by the dispersion or mixing can make it 
greater than the kinematic age. It is clear that apparent age increases as the 
time for the chemical measure increases. 



In a hypothetical one-dimensional flow field, groundwater velocity (v) is as- 
sumed to be 100 m/yr with the dispersion coefficient of 1000m. For the sake 
of simplicity, activity at time zero is assumed to be Gaussian with the peak 
Sit X = (to — lyr). Fig. shows the spatial distribution of tritium activ- 
ity at i = 0, 10, and 20 years and Fig. [6]3 shows the yearly time series at 
X = 1000m. The spatial distribution of tritium activity a.t t = lOyears indicates 
that it peaks at a; = 1000m and thus the kinematic age or the mean age can 
be calculated as 10 years. This is valid also for the conservative tracer. If the 
simple rate of decay is used to calculate the apparent age at t = lOyears, then 
Aa{x = Xpeak) ~ 31.5. Usiug the ratio between non-decaying and decaying 
solutes, one obtains Aa{x = Xpeak) — — (1/A) In (C/Ca=o) ~ ll.Oyears. 



Fritz et al. ( 1991 1 estimated groundwater mean residence time using tritium 
activities measured from the samples collected at the same location in different 
years. They suggested to use the criterion p{t) = C{t) exp{—Xt)/C{t + At) to 
determine whether the peak passed the location or not. If p{t) > 1, then the 
peak has passed the observation point, whereas if p{t) < 1, then the observation 
point is sampling the leading edge. From Fig. [6]d, p{t « 8.8 years) = 1 and so 
the age is determined to be 8.8 years at x = 1000 m. The method by [Fritz et al. 



( 1991 1 is derived from the decay equation. If the observation point is sampling 
the leading edge, then dC{t)/dt > ~XC{t) and for the tail, dC{t)/dt < ~\C{t) 
since the advection and dispersion increases the concentration before the peak 
passes the location while concentration decreases after the peak. Although this 
is a semi-quantitative approach to determine the relative location of the peak in 
the flow system, it closely approximates the kinematic. |Delhez and Deleersni- 



jder (2008) however show how such apparent ages systematically underestimate 



the mean age. 



It is implied from the above analyses that the isotopic age dating could pro- 



15 




Figure 6: (a): Spatial distribution of tritium activity; (b): Temporal distribution of 
tritium activity, with indication on the criterion p{t) of Fritz et al. (1991|. Solid lines: 
decaying tracer; Dashed line: conservative tracer. 



vide an approximate mean age or kinematic age when the flow system is homo- 
geneous and uniform (or approximated as one-dimensional), when the mixing is 
minimal, and if the spatial and temporal location of the peak is used to measure 
the age. However, the calculated age is based on the radioactivity (radiometric 
age) and may not be representative for the mean residence time of the sample 
or groundwater velocity at given locations. 



4. Numerical Example 



Schwartz et al. (2010) analyzed the C apparent age at Yucca Mountain in 



Nevada, in a steady flow system and also in a historic transient since the last 
glacial maximum, and indicated that accounting for the historic transients in 
recharge in the flow system makes a significant difference in apparent age distri- 
bution at present compared to that simulated using the present-day steady flow 
system. In this study, the mean age simulated using Eq. ^ is compared to the 
radiometric age translated from the results of the simulation using Eq. (U m 
the same steady flow system used by Schwartz et al. ( 2010[ ) (Fig. [t]) . Details for 
the domain, initial and boundary conditions, and all the parameters used for 
the simulations are referred to Schwartz et al. (2010): in the two-dimensional 



vertical flow domain (approximately 20 km long and 2 km deep) consisting of 
different volcanic aquifers and aquitards, water table location is 500 m below 
the ground surface or deeper, and recharge ranges between 0.05 to 10 mm/year. 



Fig. [T] shows the distributions of the simulated mean and radiometric ages in 
the steady flow system. From the results in Fig. [7] it is clear that the simulated 
ages are similar in pattern (if the water is older in terms of radiometric age, then 
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Age [years] 




Apparent Age / Mean age 




Figure 7: (a) Mean and (b) radiometric ages simulated using Eqs. ([5| and ^ in Yucca 
Mountain, Nevada, (c) Relative difTerence. Steady flow field is assumed for both 
simulations and the black solid lines indicate the location of the water table. 



the mean residence time is greater based on the mean age) but the radiometric 
age underestimates the mean residence time, especially when the water becomes 
older, agreeing with the analysis in the previous sections. In the unsaturated 
zone above the water table where transport is mostly advection-dominated, 
the two age distributions are almost identical except for the low permeable-low 



recharge region where the Caldera Barrier exists (Schwartz et al. (20101). In the 



saturated zone, the two age distributions are relatively distinct since the mixing 
of waters of different origins and the dispersion processes play more important 
roles. 



5. Conclusions 

The presented one-dimensional and two-dimensional analyses reveal that the 
differences between mean age (as obtained from the direct physical modelling 
of mean age fate) and apparent age (as obtained from isotope concentrations) 
can be very important, and that the competition between decay and dispersion 
coefficient is a major factor that can explain such differences. More precisely. 
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the larger the product XD, the larger the difference between mean and ap- 
parent ages. It is also shown that apparent age systematically overestimates 
kinematic age and underestimates mean age. This statement is verified in the 
context of aquifer-aquitard systems, for which coupled one-dimensional and two- 
dimensional solutions also show that mean age in the aquitard is significantly 
higher than the apparent age although both ages are influenced by the mixing 
processes in the aquitard. A correction to the average apparent velocity Va that 
is obtained from two apparent age dates at two distinct locations is also pro- 

2 

vided. This correction is a direct function of the dimensionless variable and 
indicates that the average velocity may be calculated from the apparent veloc- 
ity only when Va > VXD. The analytical analysis is concluded by the study of 
the effects of transverse mixing on apparent and mean ages. The results indi- 
cate a relatively minor effect of dimensionality for uniform flow flelds, and that 
apparent and mean ages are not influenced by transverse dispersion along the 
streamline from the source. 

The differences in apparent age and mean age observed with simple system con- 
flgurations can be expected to be much larger in more realistic and complex 
situations, e.g. when several flow systems tend to mix, with the consequence 
that dispersion of ages is not the result of mechanical dispersion only, but also 
a consequence of the travel times repartition within each flow system (bcnng 
itself a consequence of the shape and length of flow paths, thus of the geological 
structure of the aquifer and of the boundary conditions). A rnimerical analysis 
of age distribution at Yucca Mountain in Nevada is presented to illustrate the 
magnitude of these age differences one can expect in realistic situations. 
Otu rcsidts suggest that the use of the classical piston-flow model (and by exten- 
sion of the exponential model and others) tends to be irrelevant for age dating 
and particularly for the interpretations that are frequently carried out to derive 
average flow velocities. Age definitions based on the radioactivity of isotopes 
may not be representative of the mean residence time of the sample or for the 
groundwater velocity at given locations, and may not be suitable in many cases 
for the calibration of hydrogeological models. 
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